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Abstract 

We substantially apply the Li criterion for the Riemann hypothesis to hold. 
Based upon a series representation for the sequence {A^}, which are certain 
logarithmic derivatives of the Riemann xi function evaluated at unity, we de- 
termine new bounds for relevant Riemann zeta function sums and the sequence 
itself. We find that the Riemann hypothesis holds if certain conjectured prop- 
erties of a sequence rjj are valid. The constants rjj enter the Laurent expansion 
of the logarithmic derivative of the zeta function about s = 1 and appear to 
have remarkable characteristics. On our conjecture, not only does the Riemann 
hypothesis follow, but an inequality governing the values \n and inequalities 
for the sums of reciprocal powers of the nontrivial zeros of the zeta function. 
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Introduction 

In this paper, we reduce the verification of the Riemann hypothesis to a conjec- 
ture concerning the behavior of certain coefficients rjk which appear in the Laurent 
expansion of the logarithmic derivative of the zeta function about s = 1. Moreover, 
should the conjectured property hold, we would then derive a result stronger than the 
Riemann hypothesis itself. Namely, we would have obtained an explicit lower bound 
for a sequence {Xj} of certain logarithmic derivatives of the xi function evaluated at 
unit argument. Since the coefficients rjj can be written as particular limits involving 
the von Mangoldt function A, these quantities seem to encapsulate both number the- 
oretic and analytic information. Such information is expected to be encountered in 
any rigorous denial or verification of the Riemann hypothesis. 

We first present an overview of our approach. We shall use the Li equivalence 
for the Riemann hypothesis to hold [25]. This equivalence results from a necessary 
and sufScicnt condition that the logarithmic derivative of the function ■ClVll ~ ^)] 
be analytic in the unit disk. The function ^ is determined from the Riemann zeta 
function ( by way of the relation ^(s) = (s/2)(s — l)7r~*/^r(s/2)C(s), where F is the 
Gamma function [13, 20, 22, 38, 37]. Then the xi function satisfies the functional 
equation ^{s) — ^{1 — s). The Li equivalence states that a necessary and sufficient 
condition for the nontrivial zeros of the Riemann zeta function to lie on the critical 
line Re s = 1/2 is that a sequence of real numbers is nonnegative for every integer 
k. In this paper, we investigate the behavior of this sequence, based upon a series 
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representation previously derived [9]. Of this representation, we are able to bound 
one finite sum and make progress in bounding the other. We find that indeed the 
sequence {A^} has only nonnegative numbers, subject to the conjectured properties 
of the sequence rjj. Among other connections, the rjj can be readily related to the 
Stieltjes constants 7^. This gives one of several avenues for further investigation of the 
Tjj expansion coefficients. This paper includes a series of appendices A-L which con- 
tain various extensions of our summation estimations, occasional alternative proofs, 
series representations analogous to that for relevant to other L-functions, tabu- 
lated numerical values, derivative relations of the Riemann zeta function, alternative 
expansions of the logarithmic derivative of the Riemann zeta and xi functions, re- 
marks on the integer-order derivatives of the Dedekind zeta and xi functions, and 
other reference material. 

We stress that we do not just conjecture as to the nature of the rj/s, but provide a 
perhaps strong plausibility argument in favor. In addition, current numerical evidence 
[29] seems to fully support our conjecture. 

The sequence {A„}^i is defined by 



Then Li's criterion for the Riemann hypothesis to hold is that all {Xn}'^=i are non- 
negative [25]. We note that Li's convention for the xi function has a factor of two 
difference: ^Li(s) = 2^(s), although this is immaterial in logarithmic derivatives such 
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[s"-Mne(s)] 



(1) 
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as Xj. We also have 

^" = Ss^l<l-^)""'"«*- (2) 

and ^(0) = — C(0) = 1/2. (Hence ^Li(O) — 1.) The approximate numerical values 
for the first few A/s are: Ai ~ 0.0230957, Ag ~ 0.0923457, and A3 ~ 0.207639 (see 
Appendix D). In fact, we have Ai = —B, where S = ln2 + (1/2) Iutt — 1 — 7/2 ~ 
-0.0230957, and 7 ~ 0.5772156649 is the Euler constant. This follows from the 
logarithmic derivative [20] 

= + (3) 

us) p\s-P p) 

where p runs over all the nontrivial zeros of the zeta function. Thus ^'{s)/^{s) — 
^p{s — p)~^, which is consistent with Cuis) — Yipi^ ~ ^1 P)- general, the A^'s are 
connected to sums over the nontrivial zeros of C,{s) by way of [25] 



P L 



(4) 



By using the Laurent expansion of the zeta function about s = 1, 



1 °° (—A Y 
as) = r + E S^7n(5 - 1)", (5) 

where the 7^ are the Stieltjes constants [20, 19] , it is possible to write a closed form 
for the Aj's. The Stieltjes constants can be evaluated from the expression 



7fe — hm > — mm , 



(6) 



and several other forms have been given [19]. For instance, we have 

A2 = 1 + 7-7^ + 7r78-21n2- In TT- 271, (7) 

and 



2 + ^TT^ - 61n2 - 31n7r - I271 + 7(8 + 2(7 - 8)7 + 671] + 872 - ^C(3) 



(8) 

It is not difficult to determine that contains the term — (—1 )■'[_;/ (j — l)!]7j_i. The 

successive A^'s can be related in several different ways, including simply 



An+l — A„ H r 

nl 



(9) 

s=l 



Alternatively, in the particular case of A2, one can write A2 = 2Ai— A^+lim^^i s^" {s) / ^{s) . 
Elsewhere, we have very recently obtained a general explicit relation for A^ in terms 
of the Stieltjes constants [11]. 

Alternative representation of Li's Aj's 
Of particular importance for the purposes of this paper is an alternative formula 
for the particular sequence of logarithmic derivatives of the Riemann xi function given 
in Eq. (1). Due to the centrahty of this result, we briefly review the derivation [9] of 
the following representation. 
Theorem 1 

m=l V^J m=2 V^J ^ 
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where the constants r)j can be written as 

— ^ hm > —A(m)\n''m . (11) 



and A is the von Mangoldt function [13, 20, 22, 38, 37, 21]. Prom the expansion 
around s = 1 of the logarithmic derivative of the zeta function, 

yF4 = -i^-^)-'-T.Vp{s-ir, (12) 

we have 

InC(s) = -Hs - 1) - E ^(s - 1)^ (13) 

giving 

In ^(s) = - In 2 + In s - J In TT + In r(s/2) - ^ ^ " 1)^- (1^) 

2 p=i P 

The radius of convergence of the expansion (12) is 3, as the first singularity encoun- 
tered is the trivial zero of ({s) at s = —2. We then evaluate 



fin n-i / \ 1 



, (15) 

s=l 



using in particular the special values -0(1/2) = —7— 2 In 2 and V^^^^ (1/2) = (— l)'^+^n!(2"+^- 
l)C(n + 1) for n > 1, where ■0 — T'/F is the digamma function and V'^''^ is the 
polygamma function. Finally, the sum in Eq. (15) over j can be converted to a sum 
over m — n — j and the simple result — Z]^=i(— 1)™(^^) = 1 used, yielding Eq. (10). 

The Laurent expansion of C'/C with the form of the constants 77^ in Eq. (11) can 
be developed by applying Theorem 1 of Ref. [21]. In this case, the counting function 
A{x) — — Z^„<xA(n) = —ip{x), where ip is the Chebyshev function, and the error 



term is given by u{x) — x — ip{x). Equation (10) has been derived alternatively in 
Ref. [6] by a method connected with A. Weil's explicit formula. Our approach is 
independent, and we believe, more direct. In Appendix G we present an extension of 
Theorem 1 which accounts explicitly for the presence of the first six trivial zeros of 
the zeta function. 

Estimation of Sums 

We characterize each of the two summation terms on the right side of Eq. (10) in 
turn. Before this discussion, we emphasize that the first few A/s may be explicitly 
written, as indicated above, and directly verified to be positive. In addition, we note 
that apparently it is already known [5] that > for all n < 2.975 . . . x 10^^. 

The sum 

S,{n)^ ±{-irh]{l-2--)am), (16) 

m=2 \''V 

may be written in several equivalent ways, including 

oo n / \ -. oo 

Equation (17) results from the use of (e.g., [1]) 

oo 

(l-2-")C(n) = ^(2A;+l)-", n>2, (18) 

k=0 

and interchange of the order of the two summations in Eq. (16). It appears to be 
equally profitable to write as 
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n 



2k + 1 



1 + 



2" A;" 



(2A; + 1)^ 



:i7) 



By using integrals estimating these forms, we obtain 



Theorem 2 



5iN>-lnn+(7-l)- + -. 



(20) 



Remark. By inserting an integral representation for ( into Eq. (16), it is possible 
to obtain the sum Si in the form 



where iFi is the Kummer confluent hypergeometric function. In this equation, 
n iFi(l — n; 2; w) = L^_]^(w), where is an associated Laguerre polynomial. 

Proof of Theorem 2. We proceed on the basis of comparison to the form of Eq. 
(17), leaving the comparison to Eq. (19) to Appendix A. In addition, we relegate 
to Appendix L the application of Euler-Maclaurin summation to Si{n), which gives 
improved estimates of the linear term. 

By making a change of variable in the integral 



Si{n) = - [1+ iFi(l-n;2;i/2)-2iFi(l-n;2;i)] 
2 Jo 



dt 



1) 




(21) 



we obtain 




(22) 



Another change of variable and an integration by parts yields 



W^^/[l-(l-t)-]f4(n-l) 



(23) 



8 



Evaluation of the integration [16] then gives 



n 1 



(24) 



Since the integrand of Eq. (21) is monotonically decreasing with k, the inequahty 
(20) readily follows. 

Corolleiry With the other sum term in Eq. (10) denoted as 5'2, we have already 
obtained the inequality A„ > (n/2) Inn— (n/2)(H- In 7r + 2 ln2) + 3/2— 15*21. Similarly, 
it is possible to readily obtain an upper bound for and A„. We have, for instance, 
Si{n) < (n/2) Inn + (7 + l)n/2 — 1/2. It is possible to develop improved coefficients 
for the linear in n term, giving tighter bounds for Si{n), as shown in Appendix L. 

Remarks. (1) Since Li has very recently obtained explicit formulas for both Dirich- 
let and Hecke L-functions [26], which are analogous to Eq. (10), we may expect our 
summation estimation techniques to aid in making progress in verifying the extended 
and generalized Riemann hypotheses. In fact, we may observe that the term T^(n) 
defined in Theorem 1 of Ref. [26] can be rewritten as 



where 5*0 is defined in Eq. (A. 15) of Appendix A. The second line of Theorem 2 of 
Ref. [26] for AE(n) may be written as 



(^) = -S*! — n In 2 



if X(-1) = 1, 



— Sn — S] 



if x(-i) = -i, 



(25) 




(26) 
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where S3 is defined in Eq. (A. 18) of Appendix A. Inequality results based upon the 
use of Eqs. (25) and (26) and the estimations of Appendix A are given in Appendix 
E, while estimations of generalized Riemann zeta functions sums are presented in 
Appendix F. 

(2) In Ref. [39], the rate of growth of the sum Si{n) is conjectured. An integral 
expression for the Li constants is written and a saddle point method is applied. How- 
ever, the dominant O(nlnn) behaviour that we have demonstrated is not rigorously 
obtained. 

It is not difficult to show that the following recursion relation exists between the 
Stieltjes constants of Eq. (5) and the coefficients rjj. We have r/o = —70 = —7 and 



This equation is equivalent to the statement [C'(5)/C('5)]C('5) — C'(*) ^nd will be 
applied later on. 

It is also of interest to relate the sequence {77^} to the sequence {(7^} which occurs 
in the expansion 



The ensueing relation figures prominently in our conjecture for the behaviour of the 
coefficients rjj. One reason for the importance of the coefficients ak is their correspon- 
dence to sums of reciprocal powers of the nontrivial zeros p of the ( function [24, 31]: 




(27) 



(28) 



k=l 



^. From Eq. (4) we see that the connection between the values A„ and the 




Further discussion of the a/s is presented 
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in Appendix J. We shall next demonstrate 
Theorem 3 

Ok = (-l)Sfe-i - (1 - 'i-^)ak) + 1, k>2. 



(29) 



Proof of Theorem 3. The key to showing that Eq. (29) holds is to evaluate the 
particular successive Riemann zeta function sums 



n=k V / 



(30) 



Prom Eq. (14) and the expansion [1] 

(7-1) , (1-7), 



In 



sP 



In 2 + 



+ 



^-i)+i:^^^?t^[i-(«-i)n (31) 



n=2 



n2^ 



we have 



^ 2 2 ^ ^ 



. ^ ^ n=2 ^ 



-f^(-i)'+f:i:^(")a-y. (32) 

p=2 P n=2j=2 \J / 

Upon comparison with Eq. (28), the constant term in this equation (i.e., the coeffi- 
cient of the (1 — s)° term) is readily shown to be zero (or see Appendix B) and the 
linear term immediately gives 



(71 = — ^ + ^ + l-ln2 = Ai. 



(33) 



We then have 



_£!M^(._,,,f:tM(«)(i-.r^-E?(i-#. 

p=2 P n=2j=2 ""^ \J J 



(34) 



k=2 
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Upon reordering the two sums we obtain 

.,.(_,)v.-|KM^f:;), ,>_2. (35, 

We now need to demonstrate Eq. (30) . We first note that 

- [C(n)-l] /n-l\ - 1 ^ (n-l){n-2)---{n-k + l) 

and 

k{k -l){k-2)...{k- j)q'-^-' = y ^J;^ , |g|<l. (37) 

Equation (37) is easily verified by mathematical induction. With Eq. (36), we have 
proceeded by interchanging two summations. However, we could have equally well 
used the integral representation 

2 poo f^~^ 

and then interchanged the order of summation and integration, which is carried out 
in Appendix B. By applying Eq. (37) to Eq. (36) we obtain the succinct result 

f [C(n)-l] /n-l\ f. 1 

If we then invoke Eq. (18) we obtain successively Eqs. (30) and (29). Theorem 3 is 
proved in another manner in Appendix C. 

Since in Eq. (12) ('{s)/({s) + (s — 1)~^ is a transcendental function analytic in 
the disc |s — 1| < 3, there are an infinite number of rjj^s which are nonzero. Prom this 
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equation we have a good deal of information on sums and moments of this sequence. 

For instance, we have 

C'(2)" 



p=0 



1 + 



C(2) 



-0.430039 < 0, 



^(-l)Prjp = 1 - ln27r ~ -0.837877 < 0, 
and more generally, by taking successive derivatives, 

oo 

- 1) • • • (p - j + l)r]p = -{-lyji - 



(40a) 
(406) 



\C'{s)] 




[Cis)\ 





(41) 



s=2 



Various formulas for zeta function derivatives, and in particular concerning the loga- 
rithmic derivative occurring in this equation, are presented in Appendix H. Further- 
more, Eq. (13) gives 

oo 

^ ' (42a) 



Vp-i 



i-c(2) = -E , 

p=l p 



and 



ln2 = ^(-l)^'^. 

p=i P 

By evaluating Eq. (12) at s = 1/2 and using the functional equation we obtain 



(426) 



f (4)%p-2-5h-Hj) 



or 



p=0 



oo , Is p 



TT 



In7r-7H h31n2 

2 



(43a) 



(436) 



An equation such as (40a) shows that not all the ?7j's can be positive. In fact, we 
can argue that there are an infinite number of Ty^'s which are positive, and an infinite 
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number which are negative. This result, together with the fact that the sequence rjj 
decreases to zero with j increasing to infinity, as shown by either Eq. (40a) or (41), 
provides a Theorem 4. 

By multiplying Eq. (12) by with g > — 1 and integrating over s from to 1 we 
obtain 

oo „1 

Y^i-lfBip +l,q + l)r]p = + 1) + 7] - / 

p=o Jo 

where B is the Beta function. 

We now characterize the behaviour of the sequence {77^} as a function of its index. 
We have 
Theorem 5 

|^,-i|<k,|+2-^C(j), (45) 

and 

Conjecture 1 (i) The sequence {aj} decays with j faster than 1/3^. Then, (ii) for 
j sufficiently large, the sequence {rjj} alternates in sign and decreases in magnitude 
as approximately 1/3^. (ni) For all j > 0, the magnitudes \rjk\ satisfy |77fe| < ^2~''. 
Furthermore, (iv) the sequence {rjj} alternates in sign for all j > 0. 

Remark. Since the radius of convergence of the expansion (12) is 3, we know that 
Irjjl cannot increase faster than 1/3^ for sufficiently large j. In fact, the expression 
Tjj ~ —7 (—1/3)-' for all J > is a very good approximation. One may readily verify 
this assertion, for instance, with the sums appearing in Eqs. (40)-(43). In addition. 



+ 



s-1 



ds, (44) 
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all currently known numerical evidence [29] supports both this approximation and 
the strict sign alternation suggested in the conjecture. 

We can not currently prove all parts of Conjecture 1, but we can offer what we 
believe to be a strong plausibility argument. In the course of this argument we do 
prove Theorem 5. In addition, based upon known properties of aj, we very recently 
have proved part (iv) on the strict sign alternation of the {rjj} sequence for all values 
of J [11]. 

(i) The first nontrivial zero pi of the zeta function is known to lie on the critical 
hne and to have ordinate approximately given by 14.134725142 (e.g., [33]), i.e., pi — 
1/2 + iai. This zero, along with its complex conjugate, dominate the sum so 
that for large k we have ak ~ aj"'^, a rate of decrease much faster than 3"^^. In 
fact, depending upon whether k is even or odd, and if even, divisible or not by 4, a 
leading behaviour of (7^ is given by one of the four forms ±2q;J^/(1/4 + af)'' (for k 
even) or ±A;q;i~^/(1/4 + af)'^ (for k odd). In addition, one may argue much more 
conservatively with the expression 

--laA + - + ^Q;l = 2 ^ ( . 1 ^ x 



2 ^^V ' V2 ' ^^V 

1 j = 4m 

j = Am+l 

-1 j = 4m + 2 ^ ' 

j = 4m + 3 

by taking the approximate largest value of the binomial coefficient and ignoring the 
attenuating effect of the powers of 1/2. This will yield a form of approximately 
(7fc ~ {2/\/k){2/ai)^ , which is still a much faster decrease with k than 1/3^^. 
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For both emphasis and clarity, we restate Eq. (29) as 

77,_i = {-iy[a^ + (1 - 2-^)C(j) - 1], 3 > 2. (47) 
By appeal to Eq. (18) we have 

oo 1 1 

2-^CO-) > (1 - 2-^)C(i) - 1 = E J^^:^ > j> 2- (48) 

Applying the triangle inequality to Eq. (47) and using the left inequality in (48) gives 
Theorem 5. 

Continuing with the argument for part (ii) of conjecture 1, by part (i), (1 — 
2~"')C(j) — 1 > then dominates in the brackets in Eq. (47) and the sign alternation 
of the sequence {rjj} then follows for sufficiently large j. (iii) With the aid of the 
recursion relation (27) or by several other means, it is possible to calculate r)k for 
any desired initial set k = 1, . . . , /cq and directly verify their sign alternation and the 
stated inequality (e.g., as in Appendix D). For larger values of k, this inequality may 
hold due to the left inequaUty in (48) when combined with Eq. (47). 

Remark. It is not essential to our main purpose here, but we may comment on 
the sign pattern of the {(Tk} sequence. The initial, and in a sense typical, sign pattern 

is simply h H with cti > 0. Initially, as k takes on the respective values 

4m, 4m + 1, 4m + 2, 4m + 3, where m is a positive integer, the sign of cr^ is given by 
+, +, — , — . This pattern continues to the point where k{k — 1) > Saf. This explains 
why (746 > rather than a^e < 0. Similar considerations apply for larger values of k 
in this sequence. 
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Mainly for reference purposes, we will now indicate other possible uses of the 
recursion relation (27). It has been proved that [4] 

TT" 

which has been improved to [31] 

|,„| < Ii±(zlM, n>l. (60) 

As an illustration of the use of such results, the combination of Eqs. (27) and (49) 
gives 

Lemma 1 

|.„l < |(!l±I)l3. (-1)1 ^'^'-'^^-^r i^^-Aw. »>1. (51) 
TT" [ n (n-A;) J 

Similarly, the inequality (50) may be applied to Eq. (27), permitting, for example, 
inductive arguments on \r]n\, but we have already deduced Theorem 5. 
We are now in position to estimate the sum 

S2{n) = -j2 (""Vm-i, (52) 



m=l 



of Eq. (10). On conjecture 1, we have that the rjj^s decrease in magnitude with j and 
always alternate in sign. This means that, similar to the behaviour of the sum 5*1, 
there is a near exponential amount of cancellation in the sum 5*2. Indeed, we have 
Conjecture 2 

|^2| <37 + C2r^^/'+^ (53) 
17 



where C2 is a positive constant and e is positive and arbitrarily small. This conjecture 
is partially motivated by our discussion elsewhere [11] of the possible connection 
between the Stieltjes and Li constants and Brownian motion [5]. 

The combination of Eq. (10) and the inequalities (20) and (53) results in 
Conjecture 3 



where the approximate numerical value of the coefficient 1 + Invr + 21n2 is 3.53. 
This estimation would show, in the absence of the last two terms on the right side, 
that already for values of n exceeding 34, we would be ensured that all A^'s are 
nonnegative. In addition, as previously mentioned, for smaller values of n one has 
only to directly calculate these particular logarithmic derivatives, from either Eq. (1), 
(2), or Eq. (10) itself and verify the nonnegativeness of the Aj's (Appendix D). 



Our program for verification of the celebrated Riemann hypothesis should now be 
clear. We have invoked the Li equivalence [25], wherein it is necessary to demonstrate 
the nonnegativity of the sequence {Xn\'^=\- Our starting point has been the reformu- 
lation of the definition (1) as the series representation (10), A„ = + 5'2 + 1 — 71(7 + 
luTT-l- 21n2)/2. Some attention to the sums 5*1 and S2 [Eqs. (16) and (52)] has been 
required because they exhibit the phenomenon of exponential cancellation. That is, 
the binomial coefficient within the two summands can take on values approaching 




(54) 



Summciry and Brief Discussion 
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2"/ ^Jn. The strict sign alternation in the summands of S\ and S'2 is critical. For the 
sum 5*1 the sign alternation is explicit, while for 5*2 it has to be deduced [11], as in 
Conjecture 1. It should be mentioned that the hue of reasoning suggested here for a 
complete proof of the Riemann hypothesis does not require strict sign alternation of 
the sequence {?7j}-it is just that this result could make the estimation of the sum S'2 
much easier. In fact, very recently we have proved the strict sign alternation of this 
sequence [11]. What is next required is an argument that correlates the magnitude 
of say r]j to 77^+1. As previously mentioned, current numerical evidence [29] seems to 
support Conjecture 1 for the behaviour of the \r\j) and {o^ sequences. Certainly 
both further computation and analysis appear to be in order. 

In contrast to the Stieltjes constants 7^, whose sign pattern is not so easy to 
discern, that of the two sequences {r^fcjfclo ^^'^ {-^fcjfcli seems to be remarkably simple. 
Upon Conjecture 1 and Eqs. (29) and (48), our approach would also yield various 
inequalities for the sums (jfe of reciprocal powers of the nontrivial zeros of the zeta 
function. 

To put it very mildly, many implications could follow from the results presented. 
We may stress that the conjectured behaviour of the sequence {77^} has multiple 
implications for the von Mangoldt and Chebyshev functions, among many others. 

We may additionally stress, that should the claims of the propositions and conjec- 
tures herein indeed be valid, we would have not only verified the Riemann hypothesis 
but produced yet a stronger result. Namely, we would have developed the inequality 
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(54) for the sequence {An}. Moreover, it seems that an asymptotic version of inequal- 
ity (54) will suffice to verify the Riemann hypothesis since evidently [5] so many X/s 
are already known to be nonnegative. 

The approach of this paper suggests that verification of the Riemann hypothesis 
may be possible within analysis. Indeed, our approach may be amenable to con- 
fronting the generalized Riemann hypothesis. Some of our sum estimations carry 
over immediately to the explicit formulas for A;^(n) and AE(n) in Theorems 1 and 2 
respectively in Ref. [26]. In turn, one is left with estimating sums which contain the 
von Mangoldt function, reciprocal powers of k, powers of Ink, and Dirichlet char- 
acters. To us, this appears to be a realistic approach to the generalized Riemann 
hypothesis for Dirichlet and Hecke L-functions and the Dedekind zeta function. We 
include in Appendix E example results along this line of investigation. 
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Figure Caption 

FIG. 1. Diagram of consecutively increasing circles of convergence corresponding 
to the inclusion of the first six trivial zeros of the Riemann zeta function in the 
expansion (G.2). The location of the first nontrivial zero pi is indicated. 
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Appendix A: Estimation of the sum 5"! of Eq. (16) and of other sums 

Here we present the derivation of the inequahty (20) for the sum 5*1 of Eq. (16), 
based upon the form (19). By making a change of variable in the integral 

n-i 



we obtain 



h{n) = [ 
Jo 



— (l-- 



- 1- 



2k 



dk, 



dy 
y2- 



(Al) 



{A.2) 



We then evaluate this integral with an integration by parts, resulting in 



n 



h{n) = -<'^/>(n)+7-l + 

2 ^0 



y 



n-l 



dy' 



+ 2-", (A3) 



where ijj is the digamma function. Since the remaining integral is nonsingular and 
0(1) the inequality (20) follows. However, we may continue much further. 
By adding and subtracting 1 in the integrand in Eq. (A. 3), we may write 



n 



n 



2 2 Jo 



where 



1-y 

2 



n-l 



^ + ^[V^(n)+7]+2-", (A4) 

y ^ 



2 



n— 1 



— = -V'C^) - 7 - / 
y Ji 



^ (l-w 



\n—l 



-dw + In 2. 



(^.5) 



'1/2 W 

The integral on the right side of Eq. (A. 5) may be evaluated in multiple ways. A 
first method is to use the Gauss hypergeometric function 2F1: 



(l-w) 



n-l 



-dw 



(l-x) 



1/2 w 



Jo 1 



n-l 



+ X 



-dx 



n2 



— 2^1(1, l;n + l;-l), (A6) 
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where 



oo 



and {.)k is the Pochhammer symbol. By using a partial fractional decomposition 
of the summand in Eq. (A. 7), it is possible to show that the sum contains a term 
n2"~^ln2, leading to a possibly new reduction of the particular 2-^1 at minus unit 
argument. By applying the well known formula (e.g., [17]) 



x{x + l)---{x + N) ^f^\k J x + k' 
whose right side is simply (—l)^A^{l/x), where A is the difference operator, A/(x) = 
f{x + 1) — f{x), we obtain 



= n 



(^.9) 



fc=o \ ^ / i=o + -'■ 
Next, we proceed alternatively. By way of binomial expansion in the integrand of 

the second integral in Eq. (A. 6), we have 

/' ^\f dx = j:{-iy T + 1), (^.10) 

Jo 1 + X \ J J 



where [16] 



givmg 



P(j + 1) = i-iy ln2 + E (All) 

k=i 
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This result is in agreement with Eqs. (A. 6) and (A. 9). The sum in Eq. (A. 12) is 

easily estimated, leading to 



-1 + 



1 ^ n 
2" - 2 Ji/2 



2" -z Ji/2 w 
Using Eqs. (A.4), (A.5), (A. 12), and (A. 13) gives 



1 , 1 

^ ^ dw < 1 - — . 

- 2" 



(A13) 



n , 



h{n) > ^[V'(n)+7 + ln2-l] -1 + 2^-". 



(A14) 



The approach followed here can be used to estimate many other Riemann zeta 
function sums of interest. For instance, we have 



m=2 



■ n ' 



c(m) = E 



(A.15) 



Then one can determine 



(^.16) 



obtained by binomial expansion and term-by-term integration. On the other hand, 
by using 

/o(n)= \\ny-\^{\-yf\%, (A17) 

and an integration by parts, we have > Iq = n[ip{n) — + Since the digamma 
function satisfies [1] ip{x) —Inx — l/2x — l/12x^ -|- 0{x~^) as x — > oo, the inequality 
So{n) > n(ln n -|- 7 — 1) -|- 1 follows for n>2. 
Additionally, we may estimate 



n ' 



3=2 \JJ 



(A18) 
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Then we have the comparison integral 

"2n / 2^ 



An integration by parts gives 

h{n) = 2n f\l - (1 - 2y)"-^]^ - 2n + 1 - (-1)", (A20) 
and then we have 

73 = i-lTn (^^) - (^) + 2(-l)"[7 + ln2 + V^(n)] + tt tan |-2n+l-( 

(A21) 

One may note that the last three terms on the right side of Eq. (A. 21) enter with 
opposite signs from Eq. (26). The form of in Eq. (A. 21) may be further rewritten 
with the use of the relation [16] tlj[{l—n)/2] — ■?/'[(n+l)/2]— 7rtan(n7r/2), together with 
the doubling formula for the digamma function, ■0[(n+l)/2] = 2ip{n) —ip{n/2)—2\R 2. 
The result is 

Is{n) = 2n {[1 + (-I)IV' W - (-1)> (^) + 7 + [1 - (-1)1 In 2| - 2n + 1 - (-1)", 

{A.22) 

which we find to be very useful in Appendix E. 

Since we have developed many estimates based upon the digamma function, it 

may be useful to record another inequality for this function. By way of Binet's first 
formula [16] 

lnr(z) = (^-l/2)ln^-^+-ln27r + y^ - - + -^j --dt, (A23) 
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we have 

1 roo /I 1 \ 

iljiz)-\nz = — +/ ]e-''dt. (A.24) 

z Jo \t — IJ 

By performing manipulations on this representation, it is then possible to show that 

The right inequality in (A. 25) also follows immediately from the fact l/(e* — 1) < 
1/2 proved in Ref. [8]. 

Another way to proceed in estimating a finite alternating sum such as Si is to 
rewrite it as a contour integral [23, 15], 

(-1) = B{N + 1, -z)f{z)dz, (A26) 

where C is a positively oriented closed curve surrounding the points i,i + 1, . . . , N, 
B{x, y) — T{x)T{y) /V{x-\-y) is the Beta function, and f{z) is an analytic continuation 
of the discrete sequence f{k) to the complex plane, with no poles within the region 
surrounded by C. When the integrand decreases sufficiently rapidly toward ±ioo, the 
asymptotic evaluation of this expression can be achieved by extending the contour 
of the integral to the left and collecting the residues at the newly encountered poles. 
However, in this paper we have been interested in gaining more information than just 
an asymptotic evaluation. 
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Appendix B: Alternative evaluation of the summation (30) 

Here we perform the sum of Eq. (30) by using the integral representation of Eq. 
(38) for the zeta function. We also record two simpler sums which are useful in the 
proof of Theorem 3. 

Upon substituting Eq. (38) into the left side of Eq. (30) we find that 

C(n) fn-l\ 1 /-^ t''-^dt 



2" [k-l) 2^+\k-l)\h 



B.l] 



\k-l)\Jo sinh(V2) 
With a change of variable and evaluation of the integral [16] we have 



n=k 

Then 



(S.3) 

which is Eq. (30). 

In connection with the proof of Theorem 3, and Eq. (33) in particular, we write 

two other zeta function sums: 

g[C(»)-l]^l^^,„v^ (B.4a) 
n2» 2 2 ' ' 

n=2 ^ ^ 
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Appendix C: Alternative proof of Theorem 3 

We proceed to deduce Theorem 3 in a way similar to the proof of Theorem 5 of 
Ref. [31]. In the process, we correct a typographical error which appears in both that 
proof and Remark 5 and Eq. (1.11) of this reference. 

From the definition of the xi function in terms of the zeta function and Eq. (28), 
we have 

4 4 + il ^ 7^ = - f^-'' 

where the sum of the last two terms on the left side is given by Eq. (12) and simply 

1/s = Zljlo(l — -s)-' for |1 — s| < 1. From the expansion [16] 

oo 

= -7 + E(-l)'C(fc)(^ - lf-\ {C.2) 

k=2 

and the doubling formula satisfied by the digamma function, ip{2z) = \[ip{z) + ip{z + 
1/2)] + In 2, we obtain 

(0 = -| + p{-ir\{k + 1)(1 - 2-'=-^)(. - 1)^ - ln2. {C.3) 

The substitution of Eq. (C.3) into Eq. (C.l) and the equating of coefficients of hke 
powers of s — 1 gives again Eq. (33) for cti from the constant term and Eq. (29) from 
the rest of the terms. This gives Theorem 3, linking sums of reciprocal powers of the 
complex zeros of the zeta function with the sequence {rjj} appearing in the expansion 
(12) of the logarithmic derivative of the zeta function. 
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Appendix D: Tabulated numerical values 



k 











-0.577216 


1 


0.0230957 


0.187546 


2 


A Anno y1 CTT 

0.0923457 


-0.0516886 


6 


0.207639 


A A1 

0.0147517 


A 

4 


0.368793 


A AA /I ro A A a 

-0.00452448 


5 


A r '7 r r /I o 

0.575543 


A AA1 A A no 

0.0014468 


6 


0.827566 


C\ f\f\f\ A'~7 V A A 

-0.000471544 


7 


1.12446 


A AAA1 C CT 1 O 

0.00015518 


o 

8 


1.46576 


A AA AA CT 1 O C O 

-0.0000513452 


9 


1 or A AO 

1.85092 


A AAAA1 '^A A-[ A 

0.0000170414 


10 


2.27934 


—5.66605 X 10 


11 


o 'vr AO 

2.75036 


1 ooror ^ , iA — 

1.88585 X 10 


12 


3.26326 


a ooAcrc: ia — 7 

—6.28055 X 10 


16 


3.81724 


O AAO /I 1 ^ >■ 1A — 7 

2.09241 X 10 


14 


4.41148 


—b.972A7 X 10 


15 


r A /I TAO 

5.04508 


2.62672 X 10 


16 


5.71711 


''7 ''7 A A O A 1*^ — Q 

—7.74484 X 10 


17 


6.42658 


o roi V* -lA — Q 

2.58144 X 10 


18 


7.17248 


o r* n A A A 1 A — 10 

—8.60444 X 10 


19 


A r o '7 /I 

7.95374 


O O/^OAO - , 1A — in 

2.86808 X 10 


20 


8.76928 


A r/^AlO - , TA — 11 

—9.56012 X 10 


O "1 

21 


9.61796 


3.18668 X 10 


99 




— l.UDZZZ X lU 


23 


11.4101 


3.54072 X 10-^2 


24 


12.3513 


-1.18024 X 10-^2 


25 


13.3210 


3.93412 X 10-13 


26 


14.3179 


-1.31137 X 10-13 


27 


15.3408 


4.37124 X 10-1^ 


28 


16.3885 


-1.45708 X 10-1^ 


29 


17.4599 


4.85694 X 10-15 


30 


18.5538 


-1.61898 X 10-15 
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K 




Vk 


Q1 

oi 


19 


Do89 


5.3966 X 10 


L6 


oz 


on 

zU 


on /1 1 

8041 


1 Tnoo '7 X y in 

— 1.79887 X 10 


-16 


OO 


O 1 

zi 


ri rr OO 

958z 


IX nn/^oo 1 n — 


17 


Q A 


no 


1 om 

1301 


—1.99874 X 10 


-17 


Q 

OO 




OI OO 

3188 


D.DD247 X 10 


18 


OO 


z5 


rr>00 

DZoZ 


o oonoo in 
— z.zzUoz X lU 


-18 


QV 
O ( 


zo 


T A OO 
(4:22 


'7 /ino'7/1 1 n — 
i A\]2iA X iU 


19 


QQ 
OO 


z7 


nT /in 

9749 


— z.4o755 X iU 


-19 


on 
oy 


29 


oono 

zzOz 


o oocro^ 1A — 

8.22527 X 10 


20 


/I n 


oU 


4774 


—2.74176 X 10 


-20 




ol 


7454 


n lonin in — 


21 




o o 
00 


noo/^ 

0z36 


o c\ A p A in- 
— o.U4o4 X lU 


-21 


/I Q 

4o 


O A 

o4 


O 1 1 1 

oiii 


1 nitX/iTxy in — 


21 


/I A 
44 


OO 


d(J7z 


o oo/(oo x^ in 

—3.38488 X 10 


-22 


4t) 


OD 


m 1 o 

9113 


1 1 OOOA V * 1 A — 

1.12829 X 10 


22 


A P. 

4o 


o8 


OOO^ 

zzz7 


O T'/^AAO ^ , 1 A 

—3.76098 X 10 


-23 


/I 7 

4/ 


o9 


r /I no 

5408 


1 ocTO/^/^ xy 1 n — 

l.zoooo X iU 


23 


/I fi 
4o 


4U 


O/^ o 

8653 


/I 1 TOO T X y in 

—4.17887 X 10 


-24 


/I o 

4y 


A O 

4z 


1955 


1 OAOAC x^ 1 A — 

1.09295 X iU 


24 


bU 


o 

4o 


C O 1 1 

5311 


A £?/(010 V, 1A 

—4.64318 X 10 


-25 


Oi 


/I /I 

44 


0*71 O 

8718 


1 r /I 7 70 ^ , in — 

1.54773 X 10 


25 


OZ 


zLfi 


91 79 
Zl / z 


— o.ioyuy X lu 


-26 


OO 


47 


5671 


1.71970 X 10~ 


26 


54 


48 


9214 


-5.73232 X 10 


-27 


55 


50 


2798 


1.91077 X 10- 


27 


56 


51 


6423 


-6.36925 X 10 


-28 


57 


53 


0089 


2.12308 X 10- 


28 


58 


54 


3795 


-7.07695 X 10 


-29 


59 


55 


7542 


2.35898 X 10- 


29 


60 


57 


1331 


-7.86327 X 10 


-30 



30 



K 




Vk 


P. 1 
Oi 


oo.oioo 


O (50100 \y 1 O — 30 

z.bziuy X iu 


oz 


crn onQO 

oy.yuoy 


S 7Q(507 \y 1 0^31 

— o. 1 oby / X iu 


OO 


Di.zyoz 


O 010Q0 \y 1 O — 31 

z.yiZcsz X iu 


P. A 


DZ.Dyo4 


O 7077C; 1 o— 32 

— y./u/ ID X iu 


00 


D4.uyo / 


Q 0Qc;01 \y 1 O— 32 

o.zooyi X iu 


DO 


DO.OUoo 


1 07S(5/1 \y 1 O— 32 
— i.U ( ob4 X iU 


/ 


P.P ni (57 

oo.yiD / 


Q c;oc;/i(5 X/ 1 ri— 33 
o.Dyi34b X iU 


OO 


oo.oooi 


1 10S/10 \x 1 0~33 

— i.iyo4y X iu 


by 


RCi 7(51 Q 

by. /bio 


Q 00/IO(5 ^ 1 O— 34 

o.yy4yb x iu 




/ i.iy4z 


1 OOIRK \y 1 O— 34 

— i.ooibo X iU 


71 
/ i 


70 (5QQ7 


A /1QSS/1 X/ 1 o — 35 

4.4ooo4 X iU 


/ Z 


7/1 none; 


1 /1 70(5 1 X/ 1 n— 35 
— i.4/ybi X iU 




/O.OC)OU 


A OQOOc; NX 10^36 

4.yozuo X iu 


7/1 


7(5 007(5 

/b.yy /b 


1 (5/1 /lOO \x 1 O— 36 
— i.b44UZ X iU 


/ 1) 


7Q /l(5Q(5 

/o.4bob 


r AQr\r\K \y 1 O— 37 

o.4oUUo X iU 


7(5 


70 CiAQ A 

/ y.y4o4 


1 SO(5(5S X/ 1 o — 37 

— i.ozbbo X iU 


77 
/ / 


OI /1Q7Q 

oi.4o i 6 


(5 nssoc; X/ 1 n— 38 
b.Uooyt) X iU 


7fi 


oz.yoo ( 


O 000(5/1 XX 1 O— 38 

— z.uzyb4 X iu 


70 


OA A AOT 

c54.44o / 


(5 7(5ccc:o \x 1 O— 39 

b. / boou X iu 




Q K 0(51 7 

oo.ybi / 


o oc;t^i(5 XX 1 o— 39 
— z.zooib X iU 


oi 


Q7 /lonn 

o / .4yuu 


7 c;i 700 XX 1 o— 40 
/ .Oi / zz X iU 


oz 


so noss 
oy.Uzoo 


o c;oc;7/i XX 1 o~40 
— Z.OUO / 4 X iU 


OO 


on c;780 


8 ^^5946 X 1 n~'^^ 


84 


92.1386 


-2.78415 X 10-^1 


85 


93.7099 


9.28051 X 10-^2 


86 


95.2924 


-3.09350 X 10-^2 


87 


96.8862 


1.03117 X 10-^2 


88 


98.4912 


-3.43723 X 10-^3 


89 


100.1076 


1.14574 X 10-^3 


90 


101.7352 


-3.81911 X 10"^^ 
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K 




Vk 


yi 


1 no Q^/i 1 


i.Z / oU4 X iU 


yz 


1 nc; no/io 


— 4.Z4o4y X iU 




-LUU. UOOZj 


1 4.i4fin X in~'^^ 


94 


108.3572 


-4.71500 X 10-^^ 


95 


110.0398 


1.57166 X 10-46 


96 


111.7328 


-5.23888 X 10-47 


97 


113.4361 


1.74629 X 10-47 


98 


115.1492 


-5.82098 X 10-48 


99 


116.8719 


1.94032 X 10-48 


100 


118.6038 


-6.46775 X 10-49 



Appendix EI: Sum estimations and lower bounds pertinent to other Dirichlet functions 

Analogous to the Corollary which we have presented below Eq. (24) of the text, 
here we develop similar lower bounds appropriate for explicit formulas for Dirichlet 
and Hecke L-functions. We make substantial use of the very recent results of Ref. 
[26] , of which we need to recall some details. We relegate to the end of this first part 
of the Appendix some relations concerning elementary sums. In the second part, we 
provide independent derivations of the major results. Theorems 1 and 2, of Ref. [26]. 

Let X be a primitive Dirichlet character of modulus r, and L{s, x) the Dirichlet 
L-function of character %. The function 





where a is if x( 



1) = 1 and a is 1 ifx(— 1) = —1, satisfies the functional equation 



s,x), with a constant of absolute value one. The function 



^ (s, x) is an entire 



function of order one and has a product representation ^ (s, x) 



^(0, x) 11^(1 ~ •s/p), where the product is over all the zeros of ^{s, x)- 



We put 




n> 1. 



{E.2) 



We presume that A^(n) > for all n 



1, 2, ... if and only if all of the zeros of ^(s, x) 



are located on the critical line Re s 



1/2. Then, Li has obtained [26] 
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where 

= -E^x(^)^n-i(ln^), (^.4) 
fc=i 

-/«)=|:(;)(-W-2-K0)-||^ for X(-1) = 1. 

= E(")(-1)'2-'C(j) for X(-1) = -1, {E.5) 

and is an associated Laguerre polynomial. We recall Eq. (25) of the text, 

r^(n) = 5i -nln2 if = 1- 

= ^0-^1, if = {E.6) 

where 5*0 is defined in Eq. (A. 15) of Appendix A. Therefore, by using the summation 
estimations presented in Appendix A, we obtain 

A;, > S^{n) + - Inn + - (^In - - 1 - 21n2j + - for x(-l) = 1, 

>5;,(n) + -lnn+-(ln--lj + - for xM) = -1- (^-7) 

In accord with the discussion of the text, we conjecture that the sum 5"^ is 'small'. 
By this we mean that S-^{n) could be 0{ri) and probably even 5"-^ is 0{n^^'^~^'^), for 
e > 0. This could result from a near exponential amount of cancellation in this sum 
due to the phases present in the Dirichlet characters. 
We next introduce the function 

Us) = CBiV^/'(27r)-T + Le + ^) , (E.S) 
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where Le is the L-series associated with an eUiptic curve E over the rational numbers, 
N is the conductor, and is a constant chosen so that = 1 [7, 30]. The function 

of Eq. (E.8) is an entire function of order one and satisfies ^e{s) — w^Ei^ — s) where 
w — {—ly with r being the vanishing order of ^e(s) at s = 1/2. 
We let 



1- 1-i 



n > 1, {E.9) 



where the sum is over all zeros p of Ce{s)- All of these zeros lie on the critical line if 
and only if [26] XE{n) > for all n = 1, 2, . . .. 
Now Li [26] has obtained the explicit formula 



(E.IO) 

where 

= -i:^i'(*)ii-i(ln4). (E.ll) 

In Eq. (E.ll), = if p\N and 6(p^) = + (3^ if (p, A^) = 1, where for each 

prime number p, and (5p are the roots of the equation — OpT + p and the values 
of ap are connected with the reduction of ii^ at p [26] . 
We recall Eq. (26) of the text, so that we may write 

XE{n) = SE{n)+n (\tl ^ " 7^ +2 (l - ^ ln2) n+S^-SQ+2n+{-lf-l, {E.12) 
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where the zeta function sum ^'3 is defined in Eq. (A. 18) of Appendix A (and see 
below, Eq. (E.15)). Then, by the results of Appendix A we obtain 



ViV / 4 \ 
Xsin) > SE{n) + nlnn + nln-^+ 3+ - ln2 n- 1. {E.13) 

ZTT \ 3 / 

Again, we conjecture that the sum Se is 0(n^/^"'"'^). The values of Op include and 
±1, so that the values of b{p'^) can either be zero or include significant sign alternation 
when p\N. Similarly, for (p, A^) = 1, the roots of — a^T + p can include ±^/p and 
[±1 ± -^1 — 4pJ/2, giving various sign changes in b{p''). When E has good reduction 
at p, —2^yp < ttp < 2y/p, so that it again appears that can have significant 

changes in sign, possibly leading to much cancellation in Se- 

Concerning Eqs. (E.6) and (E.12) we record and briefiy discuss some elementary 
summation results. We have 



in2 = y ^ — -y-, 



(E.Ua) 



givmg 



leading to 



4 °° 1 



2 00 o 

-3 + i:;i(^-2(l-l"2). (^.14c) 



In addition, we have 
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= L g + E Q -(if) 

fc=3,odd fc=2,even 

= E (^) - C(J) = 2^[C(J) - 1] - C(J)- (£^-15) 
This equation is a restatement of the relation between the Riemann zeta function and 
the Hurwitz zeta function C(s, a): C{s) = ({s, 1) = (2* - l)"^C(s, 1/2)- 

The sum of Eq. (E.15) also has a close relation to ijj^'^\l/2), where ^'^^^ is the 
polygamma function, because [16] 

/I \ oo 1 

Then, with the use of Eq. (E.15), we have 

= (-l)"+in!(2'^+^ - l)C(n + 1), {E.16b) 

which is the expected result. In general, we have tjj^'^\x) = {—iy^^^n\C{n + l,x). 

We may also write an integral representation for the polygamma function which 
is very useful for evaluating terms in explicit formulas for sums over zeros of zeta 
functions. By differentiating an integral representation for ip{z) + 7, we have 

^ ^ ^ ^ Wo 1 - e-* 2 Jo sinh(i/2) ' ^ ^ 



giving the specific values 



^'"-"(i) = (-l)"'2"-7;£i^<i.. (E.lSa) 

POO qiin-lp-y 

= (-l)™2"'-iy ^———dy, {E.lSb) 
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and 
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Appendix EII: Explicit Formulas for Dirichlet and Hecke L-Functions 

Here we give alternative derivations of the very recent main results of Li [26], 
Theorems 1 and 2, of Ref. [26]. The procedure is very similar to the proof of Theorem 
1 of the text. The Riemann zeta function case extends since the Dirichlet and Hecke 
L-functions also have product expansions over their zeros and have explicit forms 
of their logarithmic derivatives. These derivations also make it very apparent that 
certain polygamma constants are the source of the elementary sums described in the 
first part of this Appendix. 

Due to the product expansion of C,{s,x), we have the formula 

= I, fc)(^ 

This equation is the analog of Eq. (15) of the text or Eq. (G.5) of Appendix G for 
the Riemann zeta function case. From Eq. (E.l) we have 

lne(5,x) = -^^1^(7) (^) +lnL(5,x), (^-20) 

giving 

where i/j — T'/T is the digamma function and A is the von Mangoldt function, such 
that A{k) = \np when A; is a power of a prime and A(A;) = otherwise. For m > 2, 
we then have 

z \ z / n 
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{E.19) 



.s=l 



where ■0^"^ is again the polygamma function. By taking the hmit s — > 1 in Eq. (E.22) 
we then obtain the representation 



A^(n) 



n 



+ E 



-E 



n 



m=2 



n\_}^_2-rn^{m-l) ( 0- + 

ymj {m — 1)! 



-ir-^^ AHx(n) , 

-^^ n=l 



(£;.23) 



„,Ji {m — i)! ~, n 
where ip{l/2) — —7 — 21n2, ip{l) — —7, 7 is the Euler constant, ip^'^~^\l) — 
(-l)'"(m - l)!C(m), and ^^"""^^1/2) is given in Eq. (E.16b). The infinite series 
in the sum S^{n) is convergent by the prime number theorem for arithmetic progres- 
sions [26, 12]. We have therefore obtained the result Eq. (E.3). 

Similarly, due to the product expansion of Ce{s), we have the formula 



A^(n) = E 



n 



1 V^J — 1)! 



(EM) 



s=l 



where from Eq. (E.8) we have 



Inisis) =lnc£ + |lnA^-sln27r + lnr(^s + i) +lnLE(^s + i) , (E.25) 



and 
ds 



1 / 1\ °° 

\n^E{s) = - IniV - ln27r + ^ (^s + -j - 



A{n)b{n) 



Re s > 1, {E.26) 



where b{n) is discussed in the first part of this Appendix. For m > 2, we then have 



^ n=l 



Hn)b{nj 



s+1/2 



In^"-^ n. 



{E.27) 
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By taking the limit s — > 1 in Eq. (E.27) we then obtain the representation 



XE{n) 



m=2 



,my (m — 1)! V2 



^ /rA (-1)"-' ~ A(n)Mn) 

whereV'(3/2) = 2(l-ln2)-7andV'^'""^H3/2) = (-l)"*(m-l)![2'"(C(m)-l)-C(m)], 
giving the result Eq. (E.12). 

Finally, we consider the case of the Dedekind zeta function (k, for which we need 
to introduce some additional notation. We let k be an algebraic number field with 
ri real places, r2 imaginary places, and degree n = ri + 2r2. The zeta function 
Cfc has the product expansion Ck{s) — 11^(1 — Np~^)~^ for Re s > 1, where the 
product is taken over all finite prime divisors of k. We put G'i(s) = Tr-'/^T{s/2) and 
^2(5) = (27r)^-"r(s), so that obviously G'i(l) = ^2(1) = 1. Then the function 

satisfies the functional equation Zk{s) = \dk\^^'^~^Zk{l — s), where dk is the discrimi- 
nant of k. 

We let Cfc = 2'^^{2'n)'"^hR/e, where h, R, and e are respectively the number of ideal 
classes of /c, the regulator of /c, and the number of roots of unity in k. With 

^k{s) = cl^s{s - l)\dkY'^Zk{s), {E.30) 

this function is entire and has ^fc(O) = 1 [40, 32]. We first present a motivation that 
an explicit formula analogous to that for and exists, and then develop the 
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corresponding explicit formula, putting 



n 



Prom Eqs. (E.29) and (E.30) we have 



In^fe(s) = -Incfc + Ins + ln(s - 1) + - In \dk\ + ri 



(£;.3i) 



— InTT — Inrf- 
2 V2 



+r2[(l - s) ln(27r) + In V{s)\ + In ^(s), 



(£;.32) 



and 



_l,j,(,) = _+_+_ta|4|+^ 



-ln7r + ^(0]+r2[-ln(27r)+^(s)]+^, Res 



where [34] 



In Eq. (E.34), p runs over the prime ideals of k and represents the norm. 
For m > 2, we then have 



(£;.33) 



(£;.34) 



(s - 1)^ 



d 



'm— 1 



m— 1 



, Re s > 1, 



(£;.35) 



where the evaluation of the first, third and fourth terms on the right side of Eq. 
(E.35) at s = 1 gives the contribution to of 



Ai'^) = E (-1)"! J {l + [(1 - 2-'")ri + r^lCM} . 



m=2 



{EM) 
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The evaluation of all of the terms on the right side of Eq. (E.33) but the second and 
last at s = 1 gives to A„ the contribution n[l + | In — |(ln7r + 7) — (ri + In 2]. 
With the aid of Eq. (E.34) we have 

°° In™ Nn 

= -EE(-ir-^^'"-^^^> Re.>l. {E.37) 

Taking the hmit s ^ 1 in Eqs. (E.33)-(E.35) and (E.37) should yield the final explicit 
representation for A„, subject to justification of the convergence of the resulting series. 



d 



m— 1 



ds 



m— 1 



Ck{s) 
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Appendix F: Further Riemann zeta function sum estimations 

As a generalization of sums such as 5*1 of Eq. (16) of the text and 5*0 of Eq. (A. 15) 
of Appendix A, we consider here sums of the form 



S,{K,n)^ {-l)"^r]n^C{m-v), (F.l) 

m=u+2 



where v + 2 < which can be extended to |Re p\ + 2 < n. In Eq. (F.l) we have 
introduced both the positive multipher k and shift v. The special cases of Eq. (F.l) 
of direct interest to this paper are i/ = with n = 1 oi n = 2^^. We are interested to 
both reformulate the sum S^, and to obtain a lower bound for it. 

If we reorder the two sums in Eq. (F.l), the inner sum takes the form 



\rnjj"^- f ^0 {n-m-u-2y.{m + u + 2y.m\ \j J 

(F.2) 

If we use the relations {m + u + 2)\ — {u + 2)\{u + 3)^ and {n — m — u — 2)\ — 
{n — u — 2)\/(2 + 1/ — ri)m: where the Pochhammer symbol {z)n — r(z + ri)/r{z), we 
obtain the terminating hypergeometric form 

n / n\ k"* / n \ 



givmg 



S.{k, n) = (-1)^ ( k'^^' ^ - 2Fi(1, 2 - n + I/; 1/ + 3; k/J). (FA) 



We now define the comparison integral 



(n \ f°° I 

v + 2rJi -p2F,{l,2-n + iy;iy + 3;K/j)dj. (F.5) 
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With the change of variable v — 1/j the integration is easily accomphshed [16] in 
terms of the generalized hypergeometric function pFg [3]: 

I,{k, n) = {-ly ^ 2) '^'^^^ 3i^2(l, l,iy + 2-n;2,iy + 3;K). (F.6) 

Equation (F.6) is simply the result of term-by-term integration and the fact that 
(l)fc/(2)fc = l/{k + 1). The relation S,^{K,n) > I^{K,,ri) then yields a family of 
inequalities. 

When K is unity we have the reduction ip{z) = — 1) 3^2(1, 1,2 — z; 2, 2; 1) — 7 
and therefore 

/.(I, n) = i-iy (^^ ^ ^ [V^(n + 1) + 7 - H,+i], {F.7) 

where ip is the digamma function and Hn — Y,k=i ^/k is the nth harmonic number 
[10], which can also be written as 

/.(I, n) = {-ly ^ ^ [H^ - H^^,]. (F.8) 

For u — 0, this gives the result of Eq. (A. 17). An alternative form of /j/(l, n) can be 
obtained by applying Theorem 1 of Ref . [36] : 



7.(l,n) = (-ir^ 



(-1/ - 1)1. 



(F.9) 



k=i k{n-iy)k 

valid for n > and integral u > —1. The relations above at unit argument can also 
be looked upon as special cases of [28] 

3^2(1, 1, i/ + 1; 2, A + 1; 1) = ^[V^(A) - V^(A - i/)], 1/ 7^ 0, Re(A - i/) > 0, 
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= AV''(A), i/ = 0, Re(A)>0, (F.IO) 

i.e., 

5I^ = ^(^)-^(^-^)' Re(A-i/)>0, A ^0,-1,-2,.... (F.ll) 

n=l '^K'^ln 

For K — 1/2 we obtain 

7.(1/2, n) = (-l)-^ + 2) 2^ 1' ^ + 2 - 2, ^ + 3; 1/2), (F.12) 

which ior 1/ — is in agreement with Iq — Ii, where Ii is given in Eq. (24) of the text. 
For K = 2 we obtain 

h{2, n) = i-iy (^^ ^ 2) =^^2(1, 1, + 2 - n; 2, 1/ + 3; 2), (F.13) 
which is in agreement with /s of Eq. (A. 22) of Appendix A when 1/ = 0. 
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Appendix G: Alternative representation of Li's A^'s 

Here we extend Theorem 1, based upon an expansion of the logarithmic derivative 
of the Riemann zeta function with a larger radius of convergence than Eq. (12) of 
the text. We demonstrate the following representation, 



m=l 



37 \5j \7j \9j UlJ \13j 



+ E (-1)™ I (1 - 2-™)CM - ^(7 + Invr + 2 ln2), (G.l) 



m=2 

where we currently do not have an arithmetic interpretation of the constants ?7j^^^. 
Prom the expansion around s = 1 of the logarithmic derivative of the zeta function, 

we have 

6 oo (12) 

InC(s) = -ln(s - 1) + Y.Hs + 2j) - E - If + constant, (G.3) 

3=1 p=l P 

giving 

In^(s) = -In2 + lns--ln7r + lnr - +constant + ^ln(s + 2j)-E 

(GA) 

With the expansion (G.2) the radius of convergence has been increased to 13 (see 
Figure 1), as we have included the contribution of all trivial zeros of ( prior to the 
encounter with the first complex zero pi. We next evaluate 

d 



I n\ 1 



i \W (^ — 1)! \_ds^' 
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(G.5) 

s=l 



using again the special values V'(l/2) = -7 -2 In 2 and ^'^"^(1/2) = (-l)"+^n!(2"+i - 
l)C(n + 1) for n > 1, where ip = F'/F is the digamma function and ip^^^ is the 
poly gamma function. Recalling the relations 

^\n{s + 2k)^- ^~P'^i;^^^- , j>l, {G.6) 

and {dp /ds^){s — l)'^ = for A; < j, and the sum 

we find Eq. (G.l). The constant term of 6 = 7 — 1 in Eq. (G.l) serves as a count 
of the number of trivial zeros of C, accounted for in the expansion (G.2) while the 
additional explicit negative terms beyond Eq. (10) of the text appearing there are 
exponentially decreasing with n. We have developed in Eq. (G.2) an expansion with 
coefficients 77]^^'' whose magnitudes increase no faster than I/IS-' for large j. 

By using Eq. (20) for Si, an extension of the Corollary of the text is 
Corollary Gl 

A„> - Inn- (l + ln7r + 21n2)- + y 



2\" /4\" /6\" /8\" /10\" /12\" 

3) + (5) + (7) + (9) + (n j + (13 



4''^ I, (G.8) 



where we have put = - Em=i (^)?7S-i- 
In developing expansions such as 

C(s- > -2k) = C'(-2A;)(s+2A;)+ic"(-2A;)(s+2A;)2+ic'"(-2A;)(s+2A;)3+0[(s+2A;)^], k>l, 

(G.Q) 
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and 

C'is) 



C{s) {s + 2k) 2C{-2k) 



it is useful to have the derivatives 



[C"(-2A;)]2 ^ C"{-2k) 



4[C(-2A;)]2 3C{-2k) 



{s + 2k) + 0[{s + 2kf], 
(CIO) 



(2n) 



C'(-2n) = (-l)-,A_Z_^(2n + 1), n> 1, 



2{2n) 



(G.U) 



and 



ln(47r^) - 2V'(2n) - - 



C(2n + l)-2C'(2n+l)|, n>l, 

(G.12) 

which follow easily by differentiating the functional equation for and putting 
s = 2n + 1. 

We note in passing the numerical value of the constant 

i ± = + 6h + ln2.) - t . -0.0926073. (G.13) 

2fc;C(-2fc) 5544 t;c(2*+i) 

We see from Eq. (G.2) that the value of 7]^ is given by 



; 2k + 



0.377918. 



(G.14) 



Therefore we can modify Eq. (G.8) to 
Corollary G2 



A„ > - Inn - (1 + ln7r + 21n2)- + — 



2\" /4\" /6\" /8\" /lO 



9/ Vll 



n / 1 O \ 'I' 



12 
13 



(12) 



n 



^ ^ ml 



m=2 



(G.15) 
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Appendix H: On derivatives of the Riemann zeta function 



Here we capture various formulas for integer order derivatives of the zeta func- 
tion. We anticipate that these could be useful in further development of the discrete 
moment problem for the coefficients rjj of Eqs. (10)-(14) and Eq. (41) and rjj of 
Eqs. (G.1)-(G.4) of Appendix G. 

We first note that the functional equation for (, along with the evaluation ({I — 
2n) = —B2n/'2n for n > 1, where are Bernoulli numbers, yields 



c'(-i; 



[H.l) 



and 

C"(-i) 



1 

6 V 8 



ln^2^ 



+ 



71 



;i-7)y + C'(2) 



In 2 



^C'(2) - ^C"(2) 



1 

12 



-l + (l-7)' + 



TT 



+ -ln2 7r + 21n7rC'(-l), 



{H.2) 



where 7 is the Euler constant, and this can be continued to higher order derivatives. 
The derived functional equation upon which Eqs. (H.l) and (H.2) are based is 

2^-T(s)C(s) cos [i^is) - ln2] + 2'-^r{s)C'{s) cos 

-|2^-T(.)C(.) sin (|.) = TT^ In 7rC(l - s) - nX'{l - s). {H.3) 

Now Elizalde [14] has given an expression for C'(— m, q), where ({z, q) is the Hurwitz 
zeta function, valid for any negative integer value of z. We state here special cases 
for q — 1: 

1 1 ^ R„, , „ 

(HA) 



1 1 °° 

c'(-i) = -^-iE 



2k+2 



6 2^^ A;(2A; + l)(2A; + 2)' 
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C'(-2) 



B 



2k+2 



36 {2k - l)k{2k + 1) {2k + 2) ' 



{H.5) 



and 

C'(-^) 



{k + iy 



-E 



-l)^(2£-l)! 

22£-1^2£ 



'min{2l-2,k) /, 

E 

r=0 



rj{2e-r-l) 



C{2i). {H.6) 



This means that we also have exphcit expressions for ('{1 + k). In particular Eqs. 
(H.l) and (H.4) or (H.6) yield an explicit form for C'(2). 

On the other hand, we can employ the integral representation, Eq. (38) of the text, 
to at least partially yield explicit values of the zeta derivatives for Re s > 1. Other 
integral representation could be used for Re s > 0, but Eq. (38) serves for illustration. 
In the following, ip denotes the digamma function and ip^^^ the polygamma function, 
as usual. 

We have 



C{s) = -^lj{s)C{s) + 



1 r°°t'-^\nt 



r(s) Jo e* - 1 



dt, 



C'{s) = -ij'{s)as)+i^\s)C{s)-2 
and therefore 



t'-^lnt 



-dt- 



1 /•°° t'-^ In^ t 



as) 



T{s) Jo e* - 1 r(s) Jo e* - 1 
~ t^-^lnt 



1 r 



-dt, 



{H.7) 



dt, {H.8) 



{H.9) 



and 



C'js) 
C{s) 



-i^'{s) + 



f'-^ln^t 



r{s)C{s)\Jo e*-l 



-dt 



T{s)C{s) 



t'-^lnt 



dt 







{H.IQ) 
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This process can be continued, 

'('{s 



mas) 



as) 

Jo e* - 1 



+ 



ns)as) 

and we note that 



■dt 



dt + 



as) " ^ ' as) 

dt / — : — —dt 



r\s)C'^{s)Jo e*-l Jo e*-l 



r^{s)e{s) 



roo t^-^\nt 

Jo e* - 1 



n 3 



-V'(^)(2) = (-iyj![C(j + i)-i]. 



(^■11) 



(i/.12) 



Following on Eq. (H.7) we have 



m=0 



m=0 



rfA'" 1 



/ 

^0 



dt. 



e* - 1 

(HAS) 



The relations of this Appendix can be developed much more for application to 
Eq. (41) of the text or elsewhere. 
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Appendix I: A digamma function integral and a Mellin transform 



Here we consider the integral 

/(A) 



(1 - 



dx 



(/.I) 



'1 x{x'^ — 1) 

of Eq. (3.4) of Ref. [6] and evaluate it in two different ways from tabulated results 
[16]. In Ref. [6] the asymptotic behaviour of this integral for large A was of interest 
for determining a certain limit denoted by the PF operation [Eq. (3.2) there]. We 
also note a polynomial of Ref. [6] that can be written as a terminating confluent 
hypergeometric series. This associated Laguerre polynomial was used in calculating 
forward and inverse Mellin transforms and we give alternative transforms. 
With the change of variable v{x) = x~^ in Eq. (I.l) we have 



Jo (1 — V 



(1 - v^) 



-dv. 



{1.2) 



Prom Ref. [16] we then obtain 



^(^ + l)+7 



A > -2, 



(7.3) 



where -0 = P'/P is the digamma function. On the other hand, we may employ a 
partial fractional decomposition in Eq. (L2) and another tabulated result [16] so that 



7(A) 



1 



2 Jo 



1 



1-^;) (l + ^;) 
where the last integral is given by [16] 



dv 



^/;(A + l)+7-ln2+ C - 

Jo 1 



+ v 



-dv 



A + 1 



(7.4) 



(7.5) 
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Use of the doubling formula V'[(A+1) /2] = 2V'(A) -'0(A/2) -2 In 2 then again yields Eq. 
(1.3). In Appendix A we have additionally given many inequalities for the digamma 
function. 

The polynomial 

was used in Ref. [6] in connection with computing Mellin transforms. By using the 



relations 

M _ (-l)'(-^).- {-n)j+, = - n)j, (7.7) 



where (.)„ is the Pochhammer symbol, this polynomial can be written as a terminating 
confluent hypergeometric series: Pn{x) — n iFi(l — n; 2; —x). In particular, a certain 
Mellin transform involving P„ converts to a Laplace transform: 

rl poo foo 

/ P„(lnx)x"-Vx= / Pn{-u)e-'''du^n iFi(l - n; 2; ii)e-""dii 

70 ^0 JO 

:F(l-n,l;2;i) = l-(l-i)". (7.8) 



n 
s 



where F is the Gauss hypergeometric function [16], which is the expected result [6]. 
In obtaining Eq. (1.8) we have used the reduction [16] 



f(l_„,l;2^_i) = (i±^. (,.9) 



Another useful point of view of the particular polynomial (1.6) is afforded by the 
theory of Laguerre polynomials L„. This family is orthogonal on the interval [0, oo) 
with decaying exponential weight function. We have the relations 

Pn{-X) ^ -^-^ ^ Ll_,{xl (/.lO) 

54 



where is an associated Laguerre polynomial. In addition, the recursion relations 
satisfied by the Laguerre polynomials [16] give 

Pni-x) = ~[Ln{x) - L^_,{x)] = -i!i±ilL,+i(a;) + ^-^^±^^L^{x). (/.II) 

Then one can recast the important Mellin transform-inverse transform pair of Eq. 
(1.8): 

r PJ-u)e-'''du^- r^^^e-'^'du^l- (l--y , Res>0. (7.12) 
^0 ^0 du \ sj 

In obtaining this equation, one can use integration by parts, the Laplace transform 
of Ln [16], and the property Ln(0) = 1. All these relations are consistent with the 
connection Ll^_-^{x) — n iFi(l — n; 2; x) and the derivative property of the confluent 
hypergeometric function. 
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Appendix J: Regeirding Conjectures 1-3 

We present here plausibility arguments in possible support of our conjectures 
concerning the detailed behaviour of the sequences {ck} and {Xj}- For this discussion, 
we let N(T) be the number of zeros of the Riemann zeta function in the critical strip 
in the upper half plane to height T. That is, N{T) denotes the number of complex 
zeros in the rectangle < Re s < 1 and < Im s < T. 

Backlund [2] showed that N{T) satisfies 



where 



T /T\ T 7 



|e(r)| < 0.137 In T + 0.443 In In T + 4.35 for T > 2. 



(J.l) 



(J.2) 



We believe then that if one were to assume certain statistical properties of the distri- 
bution of the Riemann zeros. Conjectures 2 and 3 would follow. 

Prom Eqs. (J.l) and (J.2) we can show that there is a constant Tq such that 



7r[iV(r + 1) - N{T)] < TTlnT for T > Tq. 



(J.3) 



If we write N{T) = M{T) + e(T), then 
M{T+1)-M{T) 



2n ^ V27ry 27r \2T 6T^ ^ 12T^ 



Then with Eq. (J.2) we have 



7V(T«)-iV(T)<±l„(^)+J^0.137 



21nr+ln(l + - 



1 , fT\ 1 

< — In — H . 

27r V27ry AtiT 



[JA) 



+0.866 In ln(r+l)+8.70. 
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T > 2. (J.5) 

Since ln(l + l/T) < l/T, we find 

N{T+1)- N{T) < 0.433 In T + 0.866 In ln(r + 1) + 8.407 + ( J.6) 
We also have that 

0.8661nln(r+ 1) < 0.135 In T for T > Tq, (J.7) 

8.407 < 0.431 In T for T > Tq, (J.8) 

and we may take To — 3 x 10^. These relations yield inequahty (J. 3). 

We now assume that the Riemann hypothesis holds and consider one possible 
bound that may result for the sums (7^ = J2j Pj^-, where {pj} represents the nontrivial 
zeros of the zeta function. The nontrivial zeros have the form pj = 1/2 + e+iaj, where 
bounds for e exist in the literature due to results on zero-free regions. As mentioned 
in the text, a zero pj enters the sum 0"^ along with its complex conjugate. We then 
consider the sums 



oo 



V- 1 ^ 1 [l-ln(M-l) + feln(M-l)] , w T, , ,01 

as an approximation to (7^, where [x] denotes the greatest integer contained within x. 
We have 

cc -j^ oo -j^ 

j=m"j j=[am]j<oce<j+l'^i 

~ [N{j + 1) - N{j)] ^ ~ Inj u ^>T (T^(^^ 

j=[am] ■' j=[0Cm\ •' 
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where we applied inequality (J. 3). Since the last sum in inequality (J. 10) is bounded 

by 

p Inu 1 [1 -ln([a^] - 1) + A;ln([aJ - 1)] 

we obtain inequality (J. 9). 
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Appendix K: Observations concerning a Dedekind xi function 

We note here some relations concerning a Dedekind xi function ^k, extending some 
of our earlier derivative results for the Riemann xi function [9] . 
We let k be an imaginary quadratic field of discriminant d, and 



\Cl\ F m,n=—oo 



m, n 



nt 



■.F{m, n) 



exp 



2Tlt 



F{m, n) 



, (K.l) 



where the first sum is over the inequivalent classes of positive definite integral quadratic 
forms of discriminant d. Then Li [27] has very recently shown that 

4 fOC 



4 foo 
W Jl 



{K.2) 



for all complex s, where w is the number of roots of unity. The function is entire 
and satisfies ^fc(s) = ^^(1 — s) and ^fc(O) = ^fc(l) = 2'^^hR/w, where the number of 
real places ri = 0, the number of complex places r2 = 1, = 1 is the regulator, and 
h is the number of ideal classes of k. 

The inequivalent classes of positive definite integral quadratic forms of discrim- 
inant d consist of the classes represented by the forms [18, 27] F{x,y) satisfying 
6^ — 4ac = d for either —a<b<a<coT:0<b<a — c. Now it has also 
been shown [27] that F{x,y) = ax^ + hxy + cy"^ satisfying this condition is such 
that F{m,n) > ^J\d\/2 for all integers m, n with n 7^ 0. It follows trivially that 
7rtF{m, n)lyl\d\ — 1 > 7rt/2 — 1 > for alH e [1, 00). We then have the immediate 
Proposition. For all the inequivalent classes of positive definite integral quadratic 
forms of discriminant d, we have ^jt(s) > for all real s. Furthermore, the integer- 
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order derivatives 

satisfy > ^^r all s > 1/2. The even order derivatives obey the condition 

C^k'^\s) > for all s > 1/2. Of course, as also seen by the functional equation 
for ^k{s), the odd-order derivatives vanish at s = 1/2; we have ^j^\s) ~ 

— s). This Proposition seems to mean that ^fe(s) has no zeros for real 

values of s. 

This Proposition extends some of the results of Ref. [9]. In turn, we may apply 
some of the explicit integration results obtained there in order to evaluate the function 

dt 



(KA) 



introduced in Ref. [27], where cu is the 9 series given by u!{t) — Y,^=i exp(— vrn^i). We 
have for the first term on the right side of Eq. (K.4) 

oo 1 oo -I oo -1 

^ / e'^-^-y^dt = ^ E -r(l/2, n\) = E - [l - ErfinV^)] , (K.5) 

where r{x,y) is the incomplete Gamma function [16] and r(l/2,n^7r) = — 
2n iFi(l/2; 3/2; — 7rn^)] and iFi is the confluent hypergeometric function [16, 3], 
such that iFi(l/2; 3/2; x) = {^/T:/2)Erf{-\/—x)/\/—x, where Erf is the error function 
(probability integral) [16]. 

For the second term on the right side of Eq. (K.4) we have 



E 

n=l 



\d\/2a 



-Trn^tj.— 1 



t-^dt - E 



n=l 



Ei 



2a 



- Ei{-T:v?) 



{K.6) 
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where Ei is the exponential integral [16]. We also have various elementary relations, 
including 

e-"S-Ut = nn^ jf""' e"'"- In + exp {-.rf^\l2a) In , 

(K.l) 

obtained by integration by parts, and 

l-^J\d\/2a „ 1 rTTn^^J\d\l2a 

/ e"™ nntdt=— e-"[lnM - ln(7rn^)]dt 

r-Kii^ ^f\d\l1a / I \ 

/ e^^lnWu — I exp(— Trn^y |ci|/2a) — exp(— 7rn^) J ln(7m^) 



(X.8) 



Moreover, the particular derivative values 
and 

can be evaluated in terms of infinite series with the analytic methods of Ref . [9] , and 
we note that ^^"^^(1) > for all nonnegative integers m. These special values of Eq. 
(K.IO) enter the particular logarithmic derivatives of Eq. (E.31) for A„. 
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Appendix L: Euler-Maclaurin summation applied to S'i(n) 
Here we apply Euler-Maclaurin summation to the form of the sum 5*1 (n) given in 
Eq. (17) of the text. Accordingly, we define the summand function 

f{k) = -—--l + -^——, k>0, n>2, (L.l) 

such that /(O) — n — 1 and /(oo) = 0. We can write the arbitrary integer order 
derivative of each term of Eq. (L.l). For the first term on the right side we have 

When evaluated at /c = 0, this term gives {—iy2^ j\n. In regard to the last term of 
Eq. (L.l) we have 

(^) {2klir^^~^^\2k+ir-'^' ^^-^"^ 

where is the Pochhammer symbol, and 



\dk J {n — r) 

which can be equally expressed as 

I^Aj A;" = r! (^"') A'-""' = {-lY{-n)rk''-\ {L.3c) 

Therefore we can write 

2-f4V^"(2A: + l)-" = 2" V /n (-l)^-"2^-'"(n),-^ n! 

\dk) ^ ' ^nVn) (2A; + 1)"+^-"^ (n-m)! ^ ^ 
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A version of the Euler-Maclaurin formula, given that all derivatives of / vanish at 
infinity, is 

E fin) = / f{x)dx - E -^f^-'-'KM), {L.5a) 

n=M •'^ m=l 

where Bj^ sue Bernoulli numbers, or 

oo „oo 1 ^ B 

j:f{n)= f{x)dx+-f{0)- Y: ^f^"""\^)- iL.hh) 

n=0 ■'^ ^ m=2,even 

By using Eq. (24) of the text for the integral in Eq. (L.5b) we therefore obtain 

m='2.,even 

where if) is the digamma function and 7 is the Euler constant. The sums in Eqs. (L.5)- 
(L.6) are meant in an asymptotic sense; they are highly unlikely to be convergent. 
Prom Eq. (L.6) we may obtain the successive approximations 

5i(n) = ^[V'H+7] + ^- E ^/^'"-'nO), n>l, (L.7a) 

and 

r) ry ^ R 

Si{n) = ^mn)+^] + ^- E -t/^-^-'Ho), n>3. [L.lh) 

Equation (L.7a) is expected to be a useful approximate upper bound to Si and Eq. 
(L.7b) an approximate lower bound for this sum. When n is sufficiently large, only 
the first term on the right side of Eq. (L.l) contributes in Eqs. (L.6)-(L.7), giving 
Bmf^^'~^\^)/nn\ — —'^"^~^Bjn{n/m), as m — 1 is always an odd integer. Given al- 
ternation in sign from i?2m to i?2m+2, these successive terms will change sign also. 
However, we emphasize that the sums in Eq. (L.6)-(L.7) are generally divergent. 
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